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The gravitational wave equations form a parametric resonance system during oscillatory reheating 
after inflation, when cast in terms of the electric and magnetic parts of the Weyl tensor. This is in 
direct analogy with preheating. For chaotic inflation with a quadratic potential, this analogy is exact. 
The resulting amplification of the gravitational wave spectrum begins in the broad resonance regime 
for chaotic inflation but quickly moves towards narrow resonance. This increases the amplitude and 
breaks the scale invariance of the tensor spectrum generated during inflation, an effect which may 
. . . be detectable. The parametric amplification is absent, however, in Ist-order phase transition models 

i and ' 'warm" inflationary models with continuous entropy production. 

Oh' 
(D 



Oh: 

I 



X 



I. INTRODUCTION 



Standard inflation ends with the universe in a near-vacuum state & hence reheating the universe successfully is a 
crucial requirement of inflationary scenarios. A full theory of this phoenix-like revival from the remains of the inflaton 
\l . field remains elusive given the complexities of non-equilibrium field theory on curved spacetime, although the so-called 
■ elementary theory has been known for over a decade and a half |l|,|^,|37|| . However, whereas previously the importance 
, of parametric resonance was recognised [||,^, it was not studied in the broad-resonance realm. It was essentially 
^ ■ the move into this non-perturbative regime - preheating - which lead to a radical paradigm shift demonstrating 
^\ [ that inflation may give rise to nonlinear fluctuations of sufficient magnitude to restore symmetries fl^ and offer a 
0^ way of gracefully exiting inflation jl^ . In addition, it has forced the study of backreaction phenomena p|,[T6t using 
nonperturbative methods, and the effects of scattering during the broad resonance regime |^ , ^ . 

However, the implications of oscillatory reheating (through a second order phase transition) and the resulting large 
5^ ' quantum fluctuations for local curvature (metric) fluctuations has been largely unexplored until now, and limited 
: to scalar modes p^-|30|. The evolution of the tensor (gravitational wave) spectrum has essentially been completely 
ignored, apart from a recent study of gravitational bremstrahlung generated through the interactions of the large 
quantum fluctuations in the inflaton and decay-product fields . Gravitational wave production during the bubble- 
wall collisions of a first order phase transition have been rather more studied However, the gravitational waves 
O |[ produced in these mechanisms are backreaction phenomena, since they are due to the scalar field fiuctuations, rather 
1^ ■ than the background zero mode evolution itself. 

In this paper we wish to demonstrate that there exists such an amplification of gravitational waves, essentially 
due to the oscillation of the zero mode of the infiaton during reheating. This is in addition to any gravitational 
bremstrahlung that may be produced by the associated scalar fluctuations. Further, we show that a direct analogy 
^ , exists in the treatment of preheating and gravitational wave evolution at the end of inflation if this occurs via a second 
' order phase transition, as we shall assume here. Indeed, both are governed by (approximately) Floquet systems Q In 
the case of the quadratic potential, V{(j>) = ■^m'^cjP , both preheating and gravitational wave amplification are initially 
well approximated by the Mathieu equation, hence the correspondence of the title. This "duality" is exhibited using 
the covariant Maxwell- Weyl form of the Einstein field equations (see section |l| and e.g. Q), and is partially hidden 
in the Bardeen formalism jl^ . This demonstration of gravitational wave amplification should be considered in the 
light of the controversy surrounding tensor amplification during instantaneous phase transitions. In these models, 
two different cosmologies such as de Sitter and radiation Friedmann-Lemaitre- Robertson- Walker (flrw) are matched 
across a single 3-space through appropriate junction conditions. Such discussions are indeed rather technical and 
subtle, dealing with discontinuities in the Einstein field equations, and do not model the actual physics of reheating. 
Further, there is as yet no absolute clarity as to the possibility of tensor amplification in these models, with claims 
both for and against ||ll|,|l^ amplification. 



*email: bruce@stardust.sissa.it 

'^A Floquet system is any set of linear ODE's with periodic coefficients. The solutions of such systems are characterized by 
resonance bands of exponentially growing modes, indexed by the momentum k. 
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Now, in the covariant approach to cosmology (see e.g. P, p6| , p^ ) , a fundamental vector is the four velocity of the 
fluid, which is used to perform a 3 + 1 splitting of spacetimc. The expansion, Q, is then defined via: 

e = n'^,a (1) 

where ; denotes the covariant derivative. In FLRW universes, 6 = 3S/S = 3H, where S is the scale factor and H is 
the Hubble constant. In the case of a minimally coupled inflaton field (j), we may choose our u° to be proportional to 
the timelike vector field orthogonal to the surfaces of constant 0, i.e. proportional to the covariant derivative (j):a @, 
with hab = gab + UaUb the projection tensor into the orthogonal 3-spaces. The convective time derivative is defined as 
Qa..b ^ Qa..b^^c^ arbitrary Q" "^, i.e. it is the covariant derivative projected along the four velocity. This coincides 
with the usual time derivative in flrw spacetimes. The evolution equation for the inflaton 0, with effective potential 
^eff('?^) I' is giv'^n by: 

^i + e<^ + ^^g(<^) + n0 = o (2) 

where ' = d/dcj). 11 in eq. generically represents the backreaction of quantum fluctuations on the zero mode 
evolution via a change to the effective mass of the inflaton. It may be written specifically as the polarisation operator 
p7t or given explicitly in certain approximations, such as Hartree factorisation or the large- iV limit of 0{N) 
vector models . 

The geometry of space (if the background is flrw) enters only through the expansion 0, whose evolution is given 
by the Raychaudhuri equation ||^, which in flat flrw spacetime is: 

e = -Y^' (3) 

Here k = Stt is the gravitational coupling constant. 

To illustrate preheating and hence the flrst axis of the correspondence between reheating and gravitational wave 
evolution, assume that (j) interacts with a light scalar fleld x, which itself has no self-interaction, via the Lagrangian 
interaction term ■^g^4>^x^- Consider the simplest effective potential for chaotic inflation: 

The solution of the inflaton equation of motion (H), when the frequency is an adiabatic invariant, is that of decaying 
sinusoidal oscillations, 0(f) = $(f) sin(m0t). In the absence of particle production the amplitude varies roughly as 
$ 'i^^j to the averaged expansion. The time evolution of the quantum fluctuations for each mode of the fleld 
is given by (neglecting m^) |^: 

Xfc + exfe+ (^^|^+5^$^sin2(m0t)^ Xfe = (5) 

This can be put in canonical Mathieu form if one neglects the expansion of the universe (0 = 0): 

xl! + [A(fc)-2gcos(2z)]xfe=0 (6) 



with dimensionless coefficients: 



Mk)-^ + 2,, ,, = g (7) 



and z = m^t. Depending on the size of the coupling, g, and the mass TO0, certain modes Xk will thus be amplified 
exponentially: Xk = Pfe(w0i) exp(/ife("^m0t), where pk are functions with the same period as the oscillations of the 
inflaton field and the positive ^fc*-"-* are the Floquet indices corresponding to the n-th instability band. The existence 
of resonance bands survives when the expansion of the universe is included . The parameter q divides the phase 

space into three broad classes with qualitatively different behaviour. The case g ^ 1 is well understood M and can 



^From here on we drop the e/f subscript from V{(f>)- It is implicit. 
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be treated perturbitively, the effects of expansion being important. The broad resonance case is described roughly by 
7r~^ < g < , ~ 10^, is non-perturbative and requires consideration of the backreaction of created particles on the 
zero mode (fiQ. The resonance bands are characterized by huge occupation numbers of produced particles, typically of 
order ~ The upper limit, q^, is for an expanding universe and is much lower in Minkowski spacetime (0 = 0) 
y . The wide resonance, q ^ q^, e volution is dominated by scattering effects which rapidly shut-off the exponential 



growth of the x fluctuations [po 21 



II. GRAVITATIONAL WAVES IN MAXWELL- WEYL FORM 

Gravitational waves (tensor perturbations) describe curvature fluctuations unassociated directly with matter. We 
will discuss their evolution using the elegant covariant formalism js] . 

The covariant derivative of the four velocity Ua is split into kinematic parts via: 

Ua-b = CTab + t^ab + -^Qhab - UaUb (8) 

where the shear aab = cr(ah) is traceless, LOab = ^[ab] is the vorticity and hab = Qab + Ua^b is both the projection tensor 
into, and the metric of, the hypersurfaces orthogonal to Ua in the case where uiab = ^. Ua is the acceleration and 
is caused by pressure gradients. Here indices surrounded by round brackets denote symmetrization on those indices 
while square brackets denote anti-symmetrization. 

A. The Electric and Magnetic Weyl Tensors 

A purely tensor description of gravitational waves is still partially lacking in the covariant approach | |3^ , but a 
sufhcient description was given first by Hawking (1966) in terms of the magnetic Hab and electric, Eab parts of the 
Weyl tensor, which are automatically traceless and gauge-invariant (since they vanish in exactly flrw spacetimes) 
@|, given by: 

Hac^\nab'''Cghcdu''u'' = VabcduV (9) 
Eac = CabcduV (10) 

where * is the usual dual (Hodge) operator. These definitions are completely analogous to those of the electric and 
magnetic fields in terms of the field strength. Fab, in standard electromagnetism |2^. A natural interpretation of the 
electric Weyl part is given by the geodesic deviation equations, which for the special case of a plane gravitational 
wave propagating along the direction are: 

t = E^p^P , a, (3 = 2,3 (11) 

where ^" is the connecting vector between orthonormal tetrads associated with the congruence of null geodesies ruled 
by gravitons. This means that we can directly attribute the physical effects of linear tensor perturbations, on e.g. a 
gravitational wave detector, to the electric part of the Weyl tensor. 

The evolution equations for Hab aud Eab are provided by the Bianchi identities: 

Rab[cd,e] = ^ C^'^^d = i?^!"^^'! - , (12) 

6 

with the above equivalence only holding in four dimensions. These yield the nonlinear gravitational analogues of 
the Maxwell equations |8] which take the form of two evolution equations, Eab, Hab and two constraint equations. 



divEab, divHab - see equations (Al), (A2), (A5), (A6), again just as in the case of Maxwell's equations. 

In addition we can associate a natural super-energy to gravitational waves in the covariant approach through the 
scalar: 



^This is always true for scalar fields with our choice of four velocity, since spatial gradients of a scalar commute. 
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TP jpab I rr rrab 



(13) 



which is in fact just the Ua-projected timehke component of the Bel-Robinson tensor, the super-stress-tensor for the 
gravitational field [ p5| , p6| . After expansion in eigenfunctions of the tensor Helmholtz equation, and using the relation 
= —k^/S^, the modes Hk, of the magnetic part Hat satisfy (see and appendix A for the derivation): 



Hi, 
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en,, 



^ + e + 0' + ^(M-P) 



Hk = 



(14) 



where /i and p are the relativistic energy density and pressure respectively. This is a simple decoupled equation (c.f. 
eq. H), while the shear satisfies: 



52 



9 



CTfe = 



(15) 



Finally the modes of the electric part obey an equation which is coupled to the shear (c.f. equation Al) 
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Ek 



ie(/.*+p)-Hi(A + p) 



(16) 



The shear can be eliminated from the Ek equation by another differentiation ||l3[ , but since the shear is the variable 
which appears in the covariant formula for the temperature anisotropics of the CMB, it will be convenient to retain 
it explicitly. The shear does not appear in the Eab equation when /i + p = 0, i.e. in exact de Sitter spacetime or in 
vacuum, be it Minkowski (0 = 0) or Milne (0 7^ 0). In this case there is an exact symmetry between the equations 
(^,^|) under the interchange Eab ^ Hab- This is the linearised version of the exact nonlinear rotational duality that 
exists in the full vacuum Bianchi identities p8[ | and which is the gravitational analogue of the electromagnetic duality 
in vacuum | |39[ | which lead to the Montonen-Olive conjecture and the modern progress in dualities of supersymmetry 
and string theory [Bo| . See p^,p3|] for further analysis of the gravity-electromagnetic duality in cosmology. 



1. Oscillatory dynamics in reheating 



Now let us specialise equations (p^p^ to the case of classical scalar field dynamics. Treating here only the case of 
a single scalar field, we have, using the equivalence of 4> with a perfect fluid: 



^ = K 



-4? -f y(0) 



-02 _ 



so that 



1 



(17) 



(18) 



Here ¥{(/)) is the effective potential of the scalar field. Note that if scalar perturbations become important then 
these relations will gather terms proportional to (V(/))2 |2^. This is precisely the case if one wishes to consistently 
study the production of tensor perturbations from gravitational bremstrahlung |^. However, it brings with it a host 
of complexities, since for example, eq.s (^ - ^6|) must be rederived in the presence of the backreaction of matter 
fluctuations, a highly non-trivial problem. We will not consider this issue further here. 

The important point is that with the above identifications, the equations (p^), (p^, ( p^ ) become generalizations of 
eq. (|^) . There thus exists a strong connection between the evolution governing reheating under a given potential and 
the equations governing the evolution Hk, Ek and ak particularly when V{(j)) is an even polynomial in (f>, as in chaotic 
inflation. 
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III. CHAOTIC INFLATION AND DUALITY 



Consider again the quadratic chaotic potential, Eq. (^. In addition to chaotic inflation, this describes the dynamics 
of the invisible axion, and the Polonyi and moduli fields of supergravity and string theory, with appropriate changes 
to the values of the masses, = V^" (</>)■ 

Here we will neglect the expansion of the universe as before to delineate the duality . This would of course not 
be adequate for the study of long-time oscillatory amplification, but is justified if the period of oscillations of (j) is 
small compared with typical averaged expansion times, i.e. ^ Q. However, the oscillations of the expansion rate 
should also be included in a full description ^ which will act as an additional source of resonance. This will be 
aortant in obtaining the tensor spectrum from oscillating cold dark matter (CDM) relics like the invisible axion 
The equation (Um for the modes, Hk, of the magnetic part Hab now becomes: 



+ —^$2 sin2(^^i) \Hk=0 (19) 





FIG. 1. A schematic of the instability chart for the Mathieu equation. The two diagonal lines correspond respectively to the 
zero-modes [k = 0) of the gravitational wave spectrum and the x-spectrum of eq. (4), for the case < m^. The vertical line 
near the left border of the chart is a sample spectrum for a given model (i.e. value of q). 

This is the Mathieu equation and is precisely the same as equation (^ for the evolution of the quantum fluctuations 

Xfc, with the replacement — > k-^ = 47rm^. Hence qn = tt^- The requirement that production of x bosons is more 
efficient than Hab amplification is then > iTrm^. If ^ 10^®, as required to match CMB observations then 
this is a weak constraint on the coupling g, namely g^ > 10^^^. Nevertheless, it is a constraint independent of <E> and 
hence applies to both chaotic and new inflationary models (which have quadratic potentials near the global minima) . 
If the constraint is not met, it implies that reheating occurs preferentially via production of gravitons rather than the 
^-channel ^. 

In figure 1 this situation is depicted schematically on the instability chart of the Mathieu equation: namely, the 
situation in which Qh > Ix- diagonal lines corresponding to the mode fc = 0, delimit the physical (i.e. upper) 

region of the chart. The single vertical line corresponds to a sample spectra for a given value of q. Note that for this 
value of q, the tensor spectrum has modes in the first fundamental resonance band, while there are no x niodes which 
lie in this band. Figure (2) shows a numerical integration of the spectrum as a function of time for modes in both 
stable and unstable bands. Figure (3) zooms in on the spectrum within a stable band. 



*That the expansion also oscillates can be seen from eq. (^). 

^This is not the only constraint to be met however. If the x field has moderate or strong self-interactions, the x-resonance 
is strongly suppressed |^^. This is also the case if S> 10^, due to rescattering effects In these cases, the graviton 

decay-channel could still be important and perhaps even dominant. 
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Can Hab be amplified in the broad resonance regime ? This is necessary if gravitational wave enhancement is to 
be really effective over the expansion. This requires qn > Tr~^ . Since qh = tt^^ this implies $ > 7r~^ in units of the 
Planck energy. In the case of chaotic inflation, the amplitude of oscillations goes as $ ^ 1/N, where N is the 
number of oscillations of 0, neglecting the non-equilibrium backreaction at the end of preheating which often leads to 
a sudden decrease in $ ]l6t . Thus during preheating proper, qn ir/N'^, and Hab is initially amplified in the broad 
resonance regime, but moves rapidly towards narrow resonance. 





Xk 






A 
q 


s2„,2 + 2g 

4> 


7r$2 


2 
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Consider now the equations for the electric Weyl field, i?a6, and the shear Uab- They form a partially decoupled 
linear system with time-dependent coefficients that form an approximately Floquet system. Here we ignore the decay 
of the amplitude. This is appropriate during preheating while the non-equilibrium effects and backreaction of graviton 
production are not too strong and the oscillation period of 0, (controlled by m^) is short compared with the expansion 
time scale. This gives: 

Ek+{^ + sin(m^i)') Ek^\ [$^^3 sin(2TO^t)] a,- (20) 



with 



/ ^2 3 . 5 \ 

CTfe+ f ^-K^0'+KgTO^0'jfTfe-O (21) 

The equation for the shear can again be cast in the form of the Mathicu equation with 

2 19 

Mk)a = + -^qa , qa = n—^'' . (22) 

m^A^ 19 48 

A comparison of parameters is given in table 1. 

Note that the shear always lies in a region of broader resonance than the magnetic part of the Weyl tensor because 
Qa > Qh and Acr{k) < Aff{k). Now q^ — with the requirement of broad resonance amplification of the shear 

* > ~ 0.567r-i - 0.18. 

Since it is the shear which directly determines the tensor contribution to the anisotropy of the CMB, this may allow 
one to place constraints on large-amplitude oscillatory reheating. Since power-law inflation is known to produce one 
of the strongest tensor signals during slow-roll |Q the signal to noise (due to cosmic variance and instrument) ratio 
for the tensor component in the CMB may be significantly larger than previously hoped |^]. The CMB anisotropy 
from a tensor signal is p4,p5|: 



jrj - - / Saab k'^k^dX (23) 



E 



The left hand side is a gauge- invariant measure of the anisotropy in the CMB, and is the wave vector ruling our 
past null cone. 

Now for purely tensor perturbations, Oab and Eab are related by |p^ : 

2 

C^afc = -gQCafc - Sab (24) 

we see that exponential growth of the shear implies exponential growth of the electric Weyl field, and hence the energy 
in gravitational waves, VIqw i increases exponentially, via eq. (|l3|). 
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18 ^0 



FIG. 2. The gravitational wave spectrum, Hk as a function of k and t, with the initial condition H{Q) = 0, H{0) = 10"* and 
starting in the broad resonance regime. As can be seen, the power spectrum becomes highly fc-dependent, breaking the usual 
scale-invariance of the envelope predicted by inflation. Further, the amplitude of the spectrum is exponentially enhanced by 
resonant reheating over its value during inflation. 




FIG. 3. A zoom of the stable band region 3 < A: < 5. The solution is the usual bounded oscillatory one. Note the scale on 
the z-axis (c.f. fig. 2). 

However, a counter example to the above situation occurs if there is non-thermal symmetry restoration ||l^ , [l6|l : start 
with (j) 3> Mpi, as in the chaotic inflation scenario, but with a Coleman- Weinberg type new- inflationary potential, 
which is flat near the origin. If preheating restores symmetry, a phase of new inflation begins with (/) = 0. This redshifts 
away the resonantly amplified tensor spectrum. The second stage of reheating will be much less effective in amplifying 
the tensor spectrum than the first if V{Q) ^ Ai^;- This is the generic case for non-thermal symmetry restoration and 
hence amplification will occur in the very narrow resonance region where expansion effects are expected to dominate, 
leaving an almost untouched, scale-invariant tensor spectrum. 
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IV. CONCLUSIONS 



The main resuh of this paper is that gravitational wave perturbations can be naturally amplified during a second 
order phase transition. Thus, in addition to the standard creation of a scale-invariant stochastic gravitational wave 
spectrum due to quantum fluctuations during inflation, there may be amplification of this spectrum during reheating 
which breaks the scale invariance and may significantly enhance the rms amplitude of the tensor perturbation spectrum 
with no a priori limit on the maximum wavelength affected. However the exact pattern and size of this "symmetry 
breaking" is highly model dependent, and is hence relevant to inflationary potential reconstruction attempts p2| , ^ . 
This amplification is qualitatively different from gravitational bremstrahlung since it is due to the coherent oscil- 
lation of the mean energy density and pressure of the inflaton, and is largely insensitive to the nature of the scalar 
field fluctuations. 

We have further examined chaotic inflation with a quadratic potential in detail and found a duality between 
the equations describing the growth of x fluctuations and those of the magnetic part of the Weyl tensor and the 
shear, the latter being the variable which determines the gravitational wave signal in the CMB; eq (p3|). Finally, 
the resonant enhancement of the tensor spectrum is completely missing in those models of inflation which involve 
a continuous production of entropy during inflation, such as in the "warm" inflation models [Q. Future work will 
examine the quantitative importance of expansion, quantum backreaction and scattering effects for the amplification 
of gravitational waves, and will compare the Weyl formalism with the more typical Bardeen approach. 
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APPENDIX A: DERIVATION OF THE WAVE EQUATIONS 

The Bianchi identities, eq. (p^), yield the following evolution equations psf : 

h"'ah\E'"' + J"* - 2Hj'ri"''^''P'''UbUp + h'^'a^^Eab + 

+6^;"* - 3£;,(™a*)'* - = -^(m + P)'^*'" , (Al) 

the 'S' equation, hah = gab + UaUb is the projection tensor into the hypersurfaces orthogonal to Ua, and 

h"'ah\H'"' - + 2Ea^*f]"''>''P'^UbUp + h"''a''^Hab + 

+6^"* - 3i/,,("CT*)" - Hs'^"'uj''>' = (A2) 

is the 'i? ' equation, where we have given only the perfect fiuid form, and where the covariant curl terms are defined 
as: 

J™* = ha^'^v'^^'''urH\.,d (A3) 

/"* = /?,J"r;*)'''"'u,.£:%,rf (A4) 



These are the fully nonlinear equations and despite their complexity, note the symmetries between the equations ( Al ) 
and ( [A2| ), broken only by the driving term —^{n + p)cTab in eq. (Al) which is absent in the Hab equation. In addition 



there are the gravitational 'div E' and 'div H' constraint equations: 

h' aE"-' ,dh'' s - v'^^^Uba^Hgd + 3HKuj' - i/i*V;b , (A5) 

h\H"'',dh\ + if''P''uba%Eqd - SE'sUj' - (/i + p)lu' , (A6) 
where the vorticity vector is cj" — ^rj°-^'^'^UbLdcd- 
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In linear theory about a FLRW background, to specify that we are interested in purely tensor, i.e. gravitational 
wave solutions, we require that there are no scalar or vector perturbations: 







(A7) 



Examining equations (A5, A^) and dropping all second order terms, this implies that the divergences themselves 
vanish: 



(A8) 



This ensures only purely radiative solutions, and since E°'a = H"^a = 0, by the construction of the Weyl tensor, we 
have the analogue of the transverse-traceless (TT) conditions usually imposed on tensor metric perturbations. 
The equations (Al,A2) can be linearised consistently giving the coupled set: 



and 



Eab + QEab + hl^lJb)cdeU''V''H'} + ^{fl + p)(Tab - 



Hab + QHab - hlrjb)adeU'V'E'} = 



(A9) 



(AlO) 



These can be converted to wave equations for Eab and Hab by differentiation (sec p3[). Similarly the shear, (Tab, 
satisfies a wave equation and the three variables form a partially - coupled system. Once one transforms the variables 
to momentum space via expansion in eigenfunctions, QabJ of the tensor Helmholtz equation, one obtains the ordinary 
differential equations (|l|), (|l|), (|l|). 
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